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HIGHLIGHTS 


•  SOFC  model  with  a  mixed  ionic— electronic  conductor  electrolyte  is  developed. 

•  The  modeling  results  agree  well  with  experimental  data. 

•  Effects  of  various  parameters  are  studied. 

•  The  distribution  of  oxygen  partial  pressure  in  the  MIEC  electrolyte  is  obtained. 
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A  polarization  model  for  a  solid  oxide  fuel  cell  (SOFC)  with  a  mixed  ionic  and  electronic  conductor  (MIEC) 
electrolyte  is  developed  based  on  charge  transport  equation  and  a  constant  ionic  conductivity 
assumption.  The  electrochemical  reaction  at  the  electrode  is  described  by  the  Butler-Volmer  equation 
and  the  energy  conservation  equation  is  included  as  an  additional  condition  to  complete  the  model.  The 
modeling  results  agree  well  with  experimental  data.  Utilizing  this  model  the  effects  of  key  parameters, 
including  conductivities  of  electrolyte,  electrolyte  thickness  and  cathode  exchange  current  density  on  the 
performance  of  the  SOFC  are  analyzed.  The  distribution  of  oxygen  partial  pressure  in  the  electrolyte  is 
also  obtained. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Acceptor  doped  ceria  (ADC)  is  a  promising  candidate  electrolyte 
material  for  solid  oxide  fuel  cells  (SOFC)  since  it  allows  the  SOFC  to 
operate  at  intermediate  temperature  (500-800  °C)  [1-4]  thanks  to 
its  high  ionic  conductivity  in  such  temperature  range.  However 
acceptor  doped  ceria  is  a  mixed  ionic  and  electronic  conductor 
(MIEC)  at  the  typical  SOFC  operating  conditions  and  has  consider¬ 
able  electronic  conductivity  at  low  oxygen  partial  pressure  due  to 
the  partial  reduction  of  Ce4+  to  Ce3+  [5],  resulting  in  significantly 
lower  open-circuit  potential  [6,7]. 

Mathematical  models  are  important  in  analyzing  and  studying 
the  complex  phenomena  in  fuel  cells  and  are  also  useful  tools  for 
fuel  cell  design  and  operations.  There  are  many  modeling  studies 
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on  the  SOFC  with  pure  ionic  conducting  electrolyte  (e.g.  YSZ)  [8- 
12],  but  only  limited  modeling  studies  [13-20  are  reported  for 
SOFCs  with  a  MIEC  electrolyte.  Choudhury  and  Patterson  [20]  ob¬ 
tained  expressions  for  the  current  and  voltage  of  a  SOFC  with  a 
MIEC  electrolyte  using  the  Wagner’s  theory,  and  with  these  ex¬ 
pressions,  I—V  relationship  of  the  SOFC  with  a  MIEC  electrolyte  can 
be  obtained.  Tannhauser  [19]  developed  a  model  for  a  SOFC  with  a 
MIEC  electrolyte  based  on  the  charge  transport  equation  and  the 
electric  neutrality  assumption,  and  obtained  the  I—V  relationship 
and  open-circuit  voltages.  Riess  [18]  developed  a  model  based  on 
the  charge  transport  equation  and  the  assumption  of  constant  ionic 
conductivity,  and  an  explicit  expression  of  electronic  current  was 
obtained.  However  all  of  the  above  models  were  based  on  revers¬ 
ible  electrode  and  the  electrode  reaction  kinetics  were  not  taken 
into  consideration. 

Riess,  Godickemeier  and  Gauckler  [16,17]  developed  an  elec¬ 
trode  reaction  model  for  a  SOFC  with  a  MIEC  electrolyte.  In  the 
model  the  overpotentials  at  both  electrodes  were  separated  into 
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three  parts,  i.e.  mass  diffusion  loss,  charge  transfer  loss  and 
chemical  reaction  loss.  However  several  parameters  were  used  in 
the  expressions  of  the  overpotential  and  these  parameters  were 
obtained  by  fitting  the  modeling  results  with  experiment  data. 
Wachsman  et  al.  [13,15  developed  a  continuum-level  analytical 
mode  for  a  SOFC  with  a  MIEC  electrolyte  according  to  the  transport 
equation  and  the  charge  neutrality  assumption.  The  electrode  re¬ 
actions  were  described  by  combining  the  Butler-Volmer  equation 
and  the  defect  equilibrium  equation.  In  this  model,  only  one 
parameter  I<  (the  equilibrium  constant)  need  to  be  calculated  by 
comparing  numerical  results  with  experimental  data. 

In  this  work  a  polarization  model  of  a  SOFC  with  a  MIEC  elec¬ 
trolyte  is  developed  according  to  the  transport  equation  and  the 
constant  ionic  conductivity  assumption,  and  the  electrochemical 
reaction  at  the  electrode  is  described  by  the  Butler-Volmer  equa¬ 
tion.  The  energy  conservation  equation  as  a  supplementary  equa¬ 
tion  is  first  introduced  to  complete  the  model.  No  model 
parameters  need  to  be  obtained  by  fitting  experiment  data.  Thus 
this  model  is  more  rigorous  than  the  existing  models  and  has  wider 
applications.  The  modeling  results  are  compared  with  experi¬ 
mental  data  and  the  effects  of  some  key  parameters  on  the  per¬ 
formance  of  SOFC  with  a  MIEC  electrolyte  are  studied. 

2.  Model  development 

A  ID  configuration  of  a  SOFC  with  a  MIEC  electrolyte  is 
considered  and  the  charge  transport  is  shown  in  Fig.  1.  In  this  model 
oxygen  ions  and  electrons  are  the  only  mobile  charges,  and  the 
electrochemical  reaction  is  assumed  to  occur  only  in  the  interface 
between  electrodes  and  electrolyte.  The  current  direction  from  the 
anode  to  cathode  is  defined  as  positive,  thus  the  ionic  current  is 
positive  and  electronic  current  is  negative.  This  opposite  direction 
is  due  to  the  different  driving  force  for  ions  and  electrons.  The 
potential  at  the  anode  is  specified  to  be  zero  and  that  at  the  cathode 
is  specified  to  be  the  cell  voltage,  V. 


2.1.  Charge  transport  in  a  MIEC 


The  fundamental  equations  of  the  model  are  the  charge  trans¬ 
port  equation  (Nernst-Planck  equation)  and  the  continuity  equa¬ 
tion  at  steady  state, 


m  =  -d, 


dcKx) 

dx 


giW  dtp (X) 

Zj-e  6x 


(1) 


a  0 


(2) 


where  J;(x)  is  the  mole  flux  of  charge  i,  D  is  the  diffusivity,  c  is  the 
charge  concentration,  a  is  the  conductivity,  zx  is  the  charge  number 
and  0  is  the  static  potential.  As  in  some  existing  studies  [21-23],  the 
assumption  of  constant  ionic  conductivity  is  used,  thus  the  con¬ 
centration  of  oxygen  ions  is  constant  in  the  electrolyte, 


(3) 


Then  the  ionic  current  density  is  obtained  from  Eqs.  (1)  and  (3) 
as: 


90 

J°2-  =  “V 


(4) 


Since  the  electronic  conductivity  of  a  MIEC  obeys  the  -1/n  po¬ 
wer  of  the  oxygen  partial  pressure  p0l  [24],  and  n  is  assumed  to  be  4 
for  acceptor  doped  ceria,  thus: 


=  ^ef>o21/4  (5) 

where  a °  is  a  temperature-dependent  coefficient  that  is  indepen¬ 
dent  of  oxygen  partial  pressure.  The  relationship  between  the  ionic 
and  electronic  current  densities  is  obtained  by  an  early  theoretical 
model  [25]  as: 


Jo2-  je  _  RTd (In  po2) 

aQ2-  ae  ~  4 F  dx 


where  R  is  the  gas  constant,  T  is  the  absolute  temperature  and  F  is 
the  Faraday  constant. 

Integrating  Eq.  (6)  and  the  relationship  between  the  ionic  and 
electronic  current  densities  can  be  obtained  as: 


je  = 


h2- 

<702-  Mq  —  1 


[M0(P°o2rV4 


(7) 


where  p^2  and  pl0i  are  the  oxygen  partial  pressures  in  the  electro¬ 
lyte  at  x  =  0  and  x  =  L,  respectively,  and 


M0  =  exp 


Lk 2_l 

RTa02- 


By  defining  equivalent  electronic  conductivity  ae  as, 


ae  = 


a 


o 


M0-l 


and  substituting  Eq.  (4)  into  Eq.  (7),  the  electronic  current  density 
can  be  expressed  as: 


^  90 
Je  ~  ~ffedx 

From  Eq.  (6),  the  distribution  of  the 
the  electrolyte  can  be  obtained  as  [25], 

[Po2(*)]-1/4  =  [(pS2)-1/4-A0]exp^ 


(8) 


oxygen  partial  pressure  in 


F  J o2~ 
RT  aQ2- 


+  A) 


(9) 


where 


Fig.  1.  The  schematic  of  the  charge  transport  in  the  SOFC  with  a  MIEC  electrolyte. 
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2.2.  The  electrochemical  reaction  in  the  electrodes 


The  electrochemical  reaction  in  the  cathode  is  described  by  the 
Butler-Volmer  equation, 


nF 

Jo2-  =  Jclc^  exp  ^  Vc )  ~  exp 


M  \nF 

-(1  -“c)^c 


Pk 

Phi 


(10) 


where  jo_c  is  the  exchange  current  density  for  the  cathode,  ac  is  the 
cathode  transfer  coefficient  (ac  =  0.5),  is  the  activation  polari¬ 
zation  at  the  cathode,  is  the  oxygen  partial  pressure  in  the 
cathode  and  p^  is  the  reference  oxygen  partial  pressure.  The 
activation  polarization  is  defined  as: 


e  RT .  Pq  .  , 

4ohm  =  4fln^  +  A^' 
y02 


(17) 


According  to  the  Kirchhoff  s  law,  the  relationship  between  the 
load  current  density  and  the  ionic  and  electronic  current  densities 
is  given  as: 


ji  =  Jo2-  +J'e  (18) 

Substituting  Eqs.  (11),  (15),  (17)  and  (18)  into  Eq.  (14),  a  sup¬ 
plementary  condition  for  the  oxygen  partial  pressure  is  obtained  as: 


(19) 


Vc  =  Vrev,c-(V-0L)  (11) 

where  Vrev,c  is  the  theoretical  reversible  voltage  of  the  cathode,  and 
(j)L  is  the  static  potential  at  x  =  L. 

The  electrochemical  reaction  in  the  anode  is  assumed  to  be 
reversible  since  the  activation  polarization  in  the  anode  is  much 
smaller  than  that  in  the  cathode  at  the  intermediate  temperature. 
Thus, 


2.4.  Open-circuit  voltage  and  leakage  current  density 

The  open  circuit  is  an  important  operating  condition  for  the 
SOFC  with  a  MIEC  electrolyte.  Since  the  electronic  current  leakage 
at  open  circuit  represents  the  worst  case,  the  open-circuit  voltage 
can  be  used  to  evaluate  the  performance  of  a  MIEC  electrolyte. 
Under  open  circuit  the  load  current  density  is  zero,  thus: 


Van  ~  Po2  —  Po2  (12) 

where  pl02  is  the  oxygen  partial  pressure  at  the  anode. 

2.3.  Energy  conservation  equation 

The  voltage  output  for  the  SOFC  with  a  pure  ionic  conduction 
electrolyte  can  be  obtained  as  [26]: 


Ji  —  Jo2-  +Je  —  0-  (20) 

Substituting  Eqs.  (4)  and  (7)  into  Eq.  (20), 

v-  =  -^[m°K)‘1/4-(p^)"1/4]-  (21) 

Combining  Eqs.  (4)  and  (19),  the  relationship  of  the  open-circuit 
voltage  and  leakage  current  density  (denoted  by  the  ionic  current 
density)  is  obtained  as, 


^  —  ^Th  Pact  Pohm  Pconc  (18) 

where  Eih  is  theoretical  reversible  voltage  of  the  fuel  cell,  ^ac t  is 
activation  overpotential,  ?7ohm  is  the  ohmic  overpotential  and  ^COnc 
is  the  concentration  overpotential.  Taking  the  power  inputs,  out¬ 
puts  and  all  the  energy  losses  into  consideration  and  using  Eq.  (13), 
the  energy  conservation  equation  can  be  obtained  as: 


exp 


F  jp2~  1 
RTaQ2- 


1/4 


1  -  exp 


+  1. 


(22) 

Eqs.  (4)  and  (8)— (10)  and  (19)  above  form  the  system  of  equa¬ 
tions  for  the  polarization  mathematical  model  for  a  SOFC  with  a 
MIEC  electrolyte. 


.  c  .  / RT.  Po2\ 

Jo2-£Th  -  V-  x  (4F  npi0  J 

=h-V+jQ2-  ■  >/c  +Jo-  '  Johin  +  ^e|'»)ohnT  H4) 

where  r/^m  is  the  ohmic  loss  for  the  ion  transport  and  r]lhm  is  the 
ohmic  loss  for  the  electron  transport.  The  concentration  losses  and 
ohmic  losses  in  the  electrodes  are  neglected  in  this  model  since 
these  losses  are  very  small  under  normal  SOFC  operating 
conditions. 

According  to  Eq.  (4),  the  ohmic  loss  for  oxygen  ion  transport  in 
the  electrolyte  is  obtained  as: 

=  W-  (15) 

The  ohmic  loss  for  electron  transport  can  be  obtained  by  the 
Joule’s  law  as: 

L 

L/eKhm=Je/^dx.  (16) 

0 

By  substituting  Eqs.  (5)  and  (9)  into  Eq.  (16)  and  integrating  it, 
the  ohmic  loss  for  the  electronic  current  is  obtained  as  (the  details 
are  shown  in  Appendix), 


3.  Comparison  with  experimental  data 

The  experimental  fuel  cell  uses  samaria-doped  ceria  (SDC, 
Smo.2Ceo.803_,5)  as  the  electrolyte,  25  wt.%  NiO  +  75  wt.%  SDC  as  the 
anode  and  25  wt.%  SDC  +  65  wt.%  Smo.sSro.sCoOs^  (SSC)  as  the 
cathode. 

3.1.  Experimental  procedures 

The  Smo.2Ceo.803_,5  (SDC)  and  Smo.sSro.sCoCU-d  (SSC)  are  syn¬ 
thesized  by  the  citrate-nitrate  process.  The  composite  anode  and 
cathode  are  mixed  respectively  with  an  appropriate  amount  of 
terpineol  with  6  wt.%  ethyl  cellulose  to  form  pastes.  The  supported 
electrolyte  is  formed  by  dry-pressing  method  to  a  pellet  (20  mm  in 
diameter)  with  a  300  MPa  pressure,  and  the  anode  paste  is  spread 
on  the  pellet  to  form  a  circular  electrode  (diameter:  10  mm,  area: 
0.785  cm2).  Then  the  anode  and  electrolyte  are  co-sintered  for  4  h 
at  1200  °C.  The  cathode  paste  is  spread  on  the  other  side  of  the 
pellet  to  form  the  cathode,  and  a  silver  lead  is  adhered  to  the  fuel 
cell  pellet  by  silver  paste.  The  fuel  cell  is  sintered  for  3  h  at  1050  °C. 
The  finished  fuel  cell  pellet  is  shown  in  Fig.  2a.  During  all  the  ex¬ 
periments,  40  seem  3%  humidified  hydrogen  (controlled  by  mass 
flow  controller  Sevenstar®  CS200,  Beijing  Sevenstar  Electronics  Co. 
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Ltd.)  was  supplied  to  the  anode  and  the  cathode  is  open  to  the 
ambient  air. 

To  measure  the  conductivity  of  the  SDC,  SDC  is  dry-pressed  into 
a  pellet  with  a  diameter  of  11.4  mm  and  thickness  of  1.5  mm,  and 
sintered  for  4  h  at  1200  °C.  Silver  paste  is  spread  onto  both  sides  of 
the  SDC  pellet  as  electrodes  as  showing  in  Fig.  2b.  Two  silver  leads 
are  adhered  to  each  side  of  the  SDC  pellet  to  form  the  four-electrode 
arrangement  in  order  to  eliminate  the  lead  resistance  and  contact 
resistance.  The  conductivity  of  the  SDC  is  measured  by  the  elec¬ 
trochemical  impedance  spectroscopy  (EIS)  method  (HCP-803,  Bio¬ 
logic  SAS).  First,  the  conductivity  of  the  SDC  pallet  is  measured  in 
the  ambient  air  environment  to  obtain  the  ionic  conductivity  of  SDC 
since  its  electronic  conductivity  is  negligible  under  such  a  condi¬ 
tion.  Then  the  total  conductivity  of  the  SDC  pellet  is  measured  in  3% 
humidified  hydrogen  at  a  flow  rate  of  40  seem.  Assuming  the  ionic 
conductivity  is  the  same  in  both  cases,  the  electronic  temperature- 
dependent  coefficient  is  determined  by  Eq.  (24), 


Table  1 

Measured  SDC  conductivities  at  different  temperatures. 

Temperature/0  C  700  600  500 

Ionic  conductivity/S  cm-1  0.0332  0.0112  0.0037 

Coefficient  of  electronic  2.77E-07  1.96E-08  6.15E-10 

conductivity/S  cm  1 


<7t  -  (Jn2- 


-  ffe  =  V-  +  ffePo21/4 


(24) 


3.2.  Comparison  of  model  results  with  experimental  data 


The  SDC  conductivities  determined  by  experiment  are  listed  in 
Table  1.  All  the  values  are  the  average  of  three  repeated  tests. 

In  the  model,  the  exchange  current  density  is  determined  using 
Eqs.  (10)  and  (22)  from  the  measured  open-circuit  voltage,  the 
oxygen  partial  pressure  at  the  anode  is  obtained  by  the  local 
equilibrium  of  the  water  vapor, 


Pu2o 
p'hj  \Jp02 


(25) 


The  modeling  results  of  polarization  curves  at  three  different 
temperatures  of  500  °C,  600  °C  and  700  °C  are  compared  with  the 
experimental  data  and  are  shown  in  Fig.  3.  The  modeling  results 
agree  reasonably  well  with  the  experimental  data  and  the  best 
agreement  is  obtained  at  600  °C. 

Fig.  4  shows  the  comparison  of  the  modeling  results  with  the 
experiment  data  for  two  electrolyte  thicknesses  at  600  °C  and 
700  °C,  respectively.  The  comparison  shows  that  the  modeling  re¬ 
sults  agree  well  with  the  experiment  data  for  both  electrolyte 
thicknesses.  Both  the  modeling  results  and  the  experimental  data 
show  that  the  performance  of  SOFCs  increases  with  the  decrease  in 
electrolyte  thickness. 


4.  Parametric  analysis  of  a  SOFC  with  a  MIEC  electrolyte 

The  characteristics  of  a  SOFC  with  a  MIEC  electrolyte  are  studied 
by  analyzing  the  effects  of  some  key  parameters,  including  the 
electrolyte  conductivities,  electrolyte  thickness  and  the  cathode 
exchange  current  density. 


4 A.  Effect  of  electrolyte  conductivity 

Fig.  5  shows  the  dependence  of  the  open-circuit  voltage  and 
leakage  current  density  on  the  ionic  conductivity  at  600  °C.  As 
shown  in  the  figure,  the  open-circuit  voltage  of  the  SOFC  with  a 
MIEC  electrolyte  slightly  increases  with  the  increase  of  the  elec¬ 
trolyte  ionic  conductivity,  and  the  leakage  current  density  increases 
approximately  linearly  with  the  ionic  conductivity.  Fig.  6  shows  the 
dependence  of  the  open-circuit  voltage  and  leakage  current  density 
on  the  electronic  temperature-dependent  coefficient  at  600  °C.  As 
shown  in  the  figure,  in  a  realistic  open-circuit  potential  range  (e.g. 
>0.75  V)  the  open-circuit  voltage  decreases  sharply  while  the 
leakage  current  density  increases  sharply  with  the  increase  of 
electrolyte  electronic  conductivity. 


4.2.  The  effect  of  electrolyte  thickness 

Fig.  7  shows  the  effects  of  electrolyte  thickness  on  the  open- 
circuit  voltage  and  leakage  current  density  at  600  °C.  The  open- 
circuit  voltage  increases  with  the  increase  of  the  electrolyte  thick¬ 
ness,  and  rate  of  increase  decreases  as  the  thickness  of  the  elec¬ 
trolyte  increases.  The  leakage  current  density  decreases  sharply 
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Fig.  3.  Comparison  of  modeling  results  with  experimental  data  at  different 
temperatures. 


Ionic  conductivity  /  S'cm'1 

Fig.  5.  The  dependence  of  the  open-circuit  voltage  and  leakage  current  density  on  the 
ionic  conductivity  at  600  °C. 


when  the  electrolyte  is  very  thin  and  the  rate  of  decrease  reduces 
very  quickly  as  the  electrolyte  thickness  increases. 

The  polarization  and  power  density  curves  for  different  elec¬ 
trolyte  thicknesses  at  600  °C  are  shown  in  Fig.  8.  With  the  increase 
of  the  electrolyte  thickness  the  performance  of  the  SOFC  exhibits 
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Fig.  6.  The  dependence  of  the  open-circuit  voltage  and  leakage  current  density  on  the 
electronic  temperature-dependent  coefficient  at  600  °C. 


significant  decrease,  and  the  performance  reduction  is  not  as  sig¬ 
nificant  when  the  electrolyte  thickness  is  greater  than  800  pm. 
However  the  performance  for  the  SOFC  with  a  thinner  electrolyte  is 
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Fig.  4.  Comparison  of  modeling  results  with  experimental  data  at  different  electrolyte  Fig.  7.  The  dependence  of  the  open  circuit  voltage  and  leakage  current  density  on  the 

thicknesses  and  different  operating  temperatures  (a.  600  °C,  b.  700  °C).  electrolyte  thickness  at  600  °C. 
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Fig.  8.  The  polarization  curves  and  power  density  curves  for  different  electrolyte  Fig.  11.  The  polarization  curves  and  power  density  curves  with  different  cathode  ex¬ 
thicknesses  at  600  °C.  change  current  densities  at  600  °C. 
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Fig.  9.  The  variations  of  the  electronic  current  density  with  the  operating  voltage  at 
600  °C. 


Fig.  12.  The  variations  of  the  electronic  current  density  with  the  operating  voltage  at 
600  °C. 
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Fig.  10.  The  dependence  of  the  open-circuit  voltage  and  leakage  current  density  on  the 
cathode  exchange  current  density  at  600  °C. 


Fig.  13.  The  distribution  of  oxygen  partial  pressure  in  the  electrolyte  under  open  cir¬ 
cuit  for  different  cathode  exchange  current  densities  at  600  °C. 
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lower  than  that  for  a  thicker  electrolyte  near  the  open-circuit 
voltage,  where  the  leakage  current  is  very  high. 

Fig.  9  shows  the  variations  of  the  electronic  current  density  with 
the  operating  cell  voltage  at  600  °C  for  different  electrolyte  thick¬ 
nesses.  The  electronic  current  density  decreases  with  the  decrease 
of  the  operating  voltage,  and  when  the  operating  voltage  is  lower 
than  a  certain  value  the  electronic  current  density  approaches  zero. 
At  the  same  operating  voltage,  the  electronic  current  density  de¬ 
creases  with  the  increase  of  the  electrolyte  thickness. 

In  summary,  the  performance  of  a  SOFC  with  a  MIEC  electrolyte 
is  significantly  affected  by  the  MIEC  thickness.  As  the  MIEC  thick¬ 
ness  increases,  the  open-circuit  voltage  increases  and  the  leakage 
current  density  decreases,  but  both  current  density  and  power 
density  at  a  given  operating  voltage  decrease.  This  is  also  true  for 
the  maximum  power  density.  Thus  the  MIEC  thickness  must  be 
optimized  to  achieve  the  desired  optimal  performance. 


current  densities  are  the  same,  which  is  consistent  with  the  results 
shown  in  Fig.  10. 

The  distributions  of  oxygen  partial  pressure  in  the  electrolyte 
under  open  circuit  for  different  cathode  exchange  current  densities 
at  600  °C  are  shown  in  Fig.  13.  The  oxygen  partial  pressure  increases 
linearly  with  x,  the  distance  from  the  anode,  from  x  =  0  until  x/L  is 
greater  than  0.95,  where  it  increases  very  rapidly.  The  distribution 
of  oxygen  partial  pressure  is  mostly  independent  of  the  cathode 
exchange  current  density,  except  near  the  cathode  side  (x/L  >  0.95), 
where  the  oxygen  partial  pressure  increases  slightly  with  increase 
of  the  cathode  exchange  current  density.  On  the  cathode  surface, 
x  =  L,  the  oxygen  partial  pressure  increases  slightly  with  the  in¬ 
crease  of  cathode  exchange  current  density. 


4.4.  The  distribution  of  oxygen  partial  pressure 


4.3.  The  effect  of  the  cathode  exchange  current  density 

According  to  the  Butler-Volmer  equation  (Eq.  (10)),  the  ex¬ 
change  current  density  represents  the  electrode  performance,  thus 
the  cathode  exchange  current  density  is  taken  as  the  parameter  to 
study  the  effect  of  the  cathode  performance  on  the  output  char¬ 
acteristic  of  a  SOFC  with  a  MIEC  electrolyte.  The  modeling  results 
are  shown  in  Figs.  10-13. 

Fig.  10  shows  the  dependence  of  the  open-circuit  voltage  and 
leakage  current  density  on  the  cathode  exchange  current  density  at 
600  °C.  The  open-circuit  voltage  initially  rises  rapidly  as  the  ex¬ 
change  current  density  increases,  and  then  the  rate  of  increase 
decreases  and  curves  become  almost  flat  when  the  exchange  cur¬ 
rent  density  is  greater  than  about  200  A  m-2.  It  is  interesting  to 
note  from  Fig.  10  that  the  leakage  current  at  the  open  circuit  is 
constant,  i.e.  it  is  independent  of  the  cathode  exchange  current 
density.  This  result  appears  to  contradict  with  Eq.  (22),  which 
shows  that  the  leakage  current  density  depends  on  the  open-circuit 
voltage.  The  reason  for  this  seemingly  contradiction  is  that  the  ef¬ 
fect  of  the  open-circuit  voltage  on  the  leakage  current  density  is 
very  slight  under  typical  SOFC  operating  condition  (0.6  <  V0c  <  hi 
and  T  <  1073),  since, 


The  distributions  of  oxygen  partial  pressure  in  the  electrolyte  for 
different  operating  voltages  at  600  °C  are  shown  in  Fig.  14.  Same  as 
at  open  circuit,  the  oxygen  partial  pressure  at  all  cell  voltages  in¬ 
creases  linearly  with  the  distance  from  anode  side,  and  when  x/L  is 
greater  than  about  0.95  it  increase  rapidly.  As  the  cell  operating 
voltage  decreases,  the  oxygen  partial  pressure  increases.  The  oxy¬ 
gen  partial  pressure  is  the  lowest  at  open  circuit. 


5.  Conclusion 

A  polarization  model  for  a  solid  oxide  fuel  cell  with  a  MIEC 
electrolyte  is  developed  based  on  the  charge  transport  equation 
and  the  energy  conservation  equation  is  included  as  an  additional 
condition  to  complete  the  model.  The  modeling  results  are  first 
compared  with  the  experimental  data  obtained  in  the  lab,  and  good 
agreements  are  obtained.  Then  the  model  is  used  to  study  the  ef¬ 
fects  of  various  key  parameters,  including  the  conductivities  of  the 
MIEC,  the  electrolyte  thickness  and  the  cathode  exchange  current 
density.  The  model  is  also  used  to  obtain  the  distribution  of  oxygen 
partial  pressure  in  the  electrolyte  at  different  operating  voltages. 
The  following  conclusions  may  be  obtained  from  the  modeling 
results: 


exp(-^VoC)  <  1. 

Actually  the  leakage  current  density  at  open  circuit  can  be 
expressed  as, 


Jo2- 


RT  &q2- 


In 


<P°o2 


-1/4 


+  1 


(26) 


Obviously,  the  leakage  current  density  is  dependent  on  the 
conductivities  of  the  electrolyte,  the  electrolyte  thickness  and  the 
anode  oxygen  partial  pressure. 

The  polarization  curves  and  power  density  curves  with  different 
cathode  exchange  current  density  at  600  °C  are  shown  in  Fig.  11. 
Both  current  density  and  power  density  increase  with  the  increase 
of  the  cathode  exchange  current  density.  It  is  interesting  to  note 
from  Fig.  11  that  the  I-V  curves  shift  upward  as  the  cathode  ex¬ 
change  current  density  increases,  but  the  shape  of  the  polarization 
curves  remains  unchanged. 

Fig.  12  shows  the  results  of  the  variations  of  electronic  current 
density  with  the  cell  operating  voltage  at  different  cathode  ex¬ 
change  current  densities  at  600  °C.  At  a  specific  operating  voltage, 
the  electronic  current  density  decreases  with  the  increase  of  the 
cathode  exchange  current  density.  Flowever  the  maximum  elec¬ 
tronic  current  densities  for  all  the  different  cathode  exchange 


1)  The  leakage  current  density  increases  almost  linearly  with 
electrolyte  ionic  conductivity,  while  the  open-circuit  voltage  is 
not  significantly  affected. 


Fig.  14.  The  distribution  of  oxygen  partial  pressure  in  the  electrolyte  for  different 
operating  voltages  at  600  °C. 
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2)  The  leakage  current  density  increases  sharply  and  the  open- 
circuit  voltage  decreases  sharply  with  the  electrolyte  elec¬ 
tronic  conductivity. 

3)  The  open-circuit  voltage  increases  with  electrolyte  thickness, 
while  the  power  density  decreases  with  electrolyte  thickness. 
Thus  the  electrolyte  thickness  should  be  optimized  depending 
on  the  desired  operating  voltage  and  other  operating  conditions. 

4)  The  cathode  exchange  current  density  has  almost  no  effect  on 
the  leakage  current  density;  however  the  open-circuit  voltage 
and  fuel  cell  performance  increase  significantly  with  cathode 
exchange  current  density. 

5)  The  oxygen  partial  pressure  in  the  electrolyte  increases  with  the 
decrease  in  operating  voltage. 
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ohm  ohmic 

cone  concentration 

Abbreviation 

ADC  acceptor-doped  ceria 

EIS  electrochemical  impedance  spectroscopy 

MIEC  mixed  ionic  and  electronic  conductor 

SDC  samaria  doped  ceria 

SOFC  solid  oxide  fuel  cell 

SSC  cathode  material  Smo.5Sro.5Co03_5 

YSZ  yttria  stabilized  zirconia 

Appendix 

The  ohmic  loss  of  electron  transport  in  the  MIEC 

According  to  the  Joule’s  law  shown  in  Eq.  (16),  the  ohmic  loss  of 
electron  transport  in  the  MIEC  can  be  calculated  by: 


Nomenclature 

A  active  surface  area,  cm-2 

c\  concentration  of  species  z,  mol  m-3 

Di  diffusivity  of  species  z,  m2  s-1 
e  elementary  charge,  1.60218  x  10-19  C 

Et h  theoretical  reversible  voltage  of  fuel  cell 
F  Faraday  constant,  96485  C  mol-1 

/  current  flowing  through  the  load,  A 

Ji  mole  flux  of  the  charge  i 

ji  total  current  density  defined  by  jL  =  1 1  A,  A  cm-2 

ji  current  density  of  charge  z,  A  cm-2 

I<  equilibrium  coefficient 

L  thickness  of  the  electrolyte,  m 

Pi  partial  pressure  of  species  z,  bar 

P  power  density,  W  cm-2 

R  ideal  gas  constant,  8.314  J  (K  mol)-1 

T  absolute  temperature,  I< 

V  cell  voltage,  V 

Voc  open-circuit  voltage,  V 

Vrev,c  theoretical  reversible  voltage  in  the  cathode,  V 
x  position  variable,  m 

Zi  valence  of  charge  z, 


^ohm  *^e 


(A-l) 


Substituting  the  electronic  conductivity  by  Eq.  (5),  the  ohmic 
loss  becomes, 


„e  =  _Je_  f  dx 

°lm  «S  i  [Po2(xf 1/4 
According  to  Eq.  (9),  the  ohmic  loss  can  be  written  as, 


(A-2) 


ne  =  —— 
'/ohm 

ue 


L, 

f / 


dx 


(Po;) 1,4  «p  ( -  H  -  ■&■£*)  -  >] ' 

(A-3) 


Greek  symbol 

ac  cathode  transfer  coefficient,  0.5 

Pi  polarization  loss,  V 

a  conductivity,  S  m-1 

0  inner  electric  potential,  V 

Pi  chemical  potential  of  species  z,  J  mol-1 


Rearranging  Eq.  (A-3),  we  have, 


^ohm 


-1/ 


exP(]frfe 


Kj  dx 


(p8 


-1/4 


02) 


JeV- 

Jo  2-C°e 


(A-4) 


Superscripts 

0  the  boundary  of  electrolyte  at  the  anode 

L  the  boundary  of  electrolyte  at  the  cathode 

I  parameters  in  anode 

II  parameters  in  cathode 

ref  parameters  for  reference 

O2-  oxygen  ion 

e  electron 

Subscripts 
02  oxygen 

O2-  oxygen  ion 

e  electron 

c  cathode 

an  anode 

act  activation  of  electrode  reaction 


Integrating  Eq.  (A-4)  from  x  =  0  to  x  =  L,  the  ohmic  loss  is  ob¬ 
tained  as, 


RTt 


Vthm  =  — iTln 


(n 0  ^“1/4  jea 02-  |  jea 02-  pyn/  F  fo2-  f\ 

}P02>  j 


(P°o2)‘1/4 


(A-5) 


According  to  Eq.  (9),  the  oxygen  partial  pressure  at  x  =  L  can  be 
obtained  as, 


(Po2 


■,'4-(p8,)-'/4exp(  -^f-L 


jey- 

V  ffe 


(A-6) 
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Substituting  Eq.  (A-6)  into  Eq.  (A-5),  the  ohmic  loss  can  be 
shown  to  be, 


^ohm 


RT,  (Po2)  1/4  exP (^-^r L) 


■  — —In 


F  ■“  (P°o2rV4 


(A-7) 


Rearranging  Eq.(A-7),  we  obtain 


„e 

Pohm  —  4f  n„0  L 


Po2 


(A-8) 

From  Eq.  (4),  the  expression  of  ohmic  loss  can  be  written  as, 

(A-9) 


e  RT ,  Po2  A  , 

^hm  =  4pln^+^ 
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